



***for use with file "prq.hayes.medinavidal.1.dta"


*************TABLE 1**********

xtset stnum year

xtpmg d.gini100b d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(gini100b totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)


*************TABLE 2**********
xtpmg d.p10famincb d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(p10famincb totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)

xtpmg d.p20famincb d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(p20famincb totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)

xtpmg d.p40famincb d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(p40famincb totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)

xtpmg d.p60famincb d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(p60famincb totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)

xtpmg d.p80famincb d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(p80famincb totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)

xtpmg d.p90famincb d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(p90famincb totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)


****following Blackburne and Frank (2007)****
***DFE***
xtpmg d.gini100b d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(gini100 totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace dfe cluster(stnum)
***PMG
xtpmg d.gini100b d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(gini100 totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) full pmg
***MG***
xtpmg d.gini100b d.totwelf_pc d.cashassist_pc d.unemp_comp_pc d.incometax d.corptax d.proptax d.salestax d.deathtax d.fedrev_pc d.uniondense d.bfnom d.bfcitizen govparty_c leg_cont d.pcinc gdp elect1 d.nonwhite d.pop65, lr(gini100 totwelf_pc cashassist_pc unemp_comp_pc incometax corptax proptax salestax deathtax fedrev_pc ///
uniondense bfnom bfcitizen govparty_c leg_cont pcinc gdp elect1 nonwhite pop65) ec(ec) replace mg

hausman mg DFE, sigmamore


*******creating figures 1-2 coeff change from min to max********

***use data file "min max figures replication.dta"

graph hbar ginilr ginisr, over(var19) scheme(s1mono) ytitle(change in coefficent)

**************figure 1********

graph hbar p10lr p10sr, over(p10) scheme(s1mono) title(10th Percentile) ytitle(change in coefficent) legend(off) name(p10)

**************figure 2********

graph hbar p90lr p90sr, over(p90) scheme(s1mono) title(90th Percentile) ytitle(change in coefficent) name(p90)

graph combine p10 p90, title("") scheme(s1mono) c(1)


